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Abstract 

We argue that the Maxwehian approximation can essentially underestimate the 
rates of some nuclear reactions in hot plasma under conditions very close to ther- 
mal equilibrium. This phenomenon is demonstrated explicitly on the example of 
reactions in self-sustained DT fusion plasma with admixture of light elements X = 
Li, Be, C. A kinetic analysis shows that the reactivity enhancement results from 
non-Maxwellian knock-on perturbations of ion distributions caused by close colli- 
sions with energetic fusion products. It is found that although the fraction of the 
knock-on ions is small, these particles appreciably affect the D-l-X and T-l-X reac- 
tion rates. The phenomenon discussed is likely to have general nature and can play 
role in other laboratory and probably astrophysical plasma processes. 
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1 Introduction 



The concept of nuclear reaction rate is broadly used in high-temperature 
plasma research. This rate determines reaction yield and ultimately specific 
nuclear power released in plasma. The yield of a reaction between plasma 
species 1 and 2 is given by 

F(l+2) = anin2((Tf)i2, (1) 
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where a — 1/2 or 1 for identical or different colliding nuclei, respectively, 
nx and are species densities. The key quantity in (1) is the reaction rate 
parameter {av)i2 defined as the six-dimensional integral in velocity space 

{av)i2^ y"/"(Vl)/2''(v2)cT(|Vi - V2|)|Vi - V2|dVidV2. (2) 

Here a is the 1+2 reaction cross section, Vi and V2 are particle velocities in the 
laboratory frame, /" and f2 are unit-normalized particle velocity distribution 
functions. In a number of cases the integral (2) can be simplified. For example, 
in Maxwellian plasma (crv) 12 takes the well-known form 

where // is the reduced mass of the colliding particles, E is their kinetic energy 
in the center-of-mass frame, T is the plasma ion temperature. The Maxwellian 
approximation is a conventional tool to study plasma under conditions close 
to thermal equilibrium. At the same time, however, ion distribution functions 
in high-temperature plasma strictly speaking are not purely Maxwellian. A 
reason of non-Maxwellian deviation lies in exothermic nuclear reactions pro- 
ceeding in the plasma. These reactions generate energetic projectiles which 
during slowing-down affect the formation of ion distributions. Charged par- 
ticles slow down in the plasma mainly via peripheral (small-angle) Coulomb 
scattering by thermal ions and electrons. This mechanism does not change 
the equilibrium form of ion distribution. However, since the energy of reaction 
products can reach several MeV, close (large-angle) collisions between them 
and thermal ions can also take place in the plasma. The probability of close 
collisions is determined by amplitudes of Coulomb and nuclear scattering, 
and their interference term. Although such processes occur at rare opportu- 
nity, they can transfer in a single event a large amount of energy and produce 
fast knock-on ions. These ions increase the population of high-energy tails 
of respective distributions, so that some deviation from Maxwellian functions 
appears. Apart from energetic charged particles, reaction-produced neutrons 
can also contribute to the knock-on perturbation mechanism if the plasma is 
sufficiently dense. 

Since the Maxwellian approximation has widely been used in laboratory and 
astrophysical plasma studies, a natural question arises whether the knock- 
on perturbation of ion distribution could in some cases appreciably change 
reaction rates and affect power balance in hot plasma systems. It has been 
recognized that for conventional DT and DD fusion plasmas the answer is 
negative. The combination of three factors - sizable reaction probabilities at 
thermal energies where the majority of ion population is concentrated, small 
fractions of knock-on deuterons and tritons, and moderate energy dependences 
of the fusion cross sections in the energy range associated with these fast ions - 
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makes the above processes poorly sensitive to slight modifications of D and T 
distribution tails. This especially concerns the resonant D+T reaction whose 
cross section has a broad maximum at deep sub-barrier energies. However, in 
systems composed of nuclei with Z > 1 the situation has still been intriguing. 
In such a system strong Coulomb repulsion suppresses transmission probability 
through the potential barrier between interacting nuclei and, in the absence 
of pronoimced low-energy resonances, the behavior of reaction cross section 
becomes steep at least at sub-barrier (sub-MeV) energy range. This suggests 
that the respective reaction can be sensitive to the form of ion distribution 
tail, so that the suprathermal reaction channel induced by knock-on ions may 
become appreciable. The purpose of this letter is to investigate the possible 
enhancement of nuclear reaction rates due to knock-on perturbations of ion 
distributions. 



2 Semi-qualitative consideration 



One can reproduce such situation in a simple two-temperature model. Let 
us describe ion distributions in a plasma as superposition of two functions 
f + f. The first one is MaxweUian; it represents the behavior of bulk ions with 
density n and temperature T. The second function is introduced to model the 
ensemble of knock-on ions with density n' < n. We assume that /' also is 
MaxweUian with some temperature T' > T. Then the total reaction yield (1) 
can be presented as 

y(l+2)=n«^feX(l + A), (4) 

where Y^y^ik — Oinin2R{Ti, T2) is the thermal yield provided by bulk particles, 
while A gives the suprathermal correction caused by bulk-fast and fast-fast ion 
interactions 

niR{Ti,T2)a ris ^2) a nin2 R{Ti,T2)' ^ ' 

Here it! denotes (av) for the 1-1-2 reaction between MaxweUian species with 
different temperatures. Substituting MaxweUian distributions with tempera- 
tures 9i and 62 for /" and in (2), we find that the two-temperature rate 
parameter R{9i,92) can be reduced to the reactivity RiT^^f) for MaxweUian 
plasma (3) with some effective temperature T^jf. 

llli -f- 1112 

This allows one to easily estimate A. Let us consider, as an example, different 
D+X systems: symmetric (X = D), nearly symmetric and resonant (X = T), 
asymmetric and involving light nuclei (X = Li, Be). We assume that only the 
D distribution is distorted, while the other particles are MaxweUian. It seems 
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reasonable to set bulk temperatures of D and X nearly eqiial, Td — Tx = T, 
and neglect the contribution of the fast- fast ion interaction term in (5). Under 
these conditions (5) takes the form 



A 



riB R{T) a 



(7) 



where {3 equals 2 (X=D) or 1 (Xt^ D). Choosing T = 10 keV typical of fusion 
plasma level and varying the unknown temperature in the 50-200 keV wide 
range, we find that the ratio R{T{),T)/ R{T) changes approximately within 8- 
50 (D+D), 5xl02-5xl0-'^ (D+^Li), 10^-5x10^ (D+^Li), GxlO^-TxlO^ (D+^Be), 
8x10^-6x10^ (D+^Be). For the D+T reaction this ratio changes from 6.29 to 
6.35, i.e. proves to be nearly constant. Such invariant-like behavior with re- 
spect to T/) results from resonant nature of the D-l-T reaction, due to which its 
reactivity rapidly increases at low temperature, peaks around 60 keV and then 
becomes essentially insensitive to plasma temperature. Thus, in the D+Li and 
D-I-Be systems a very small fraction of knock-on deuterons n^/nD < 0.1% 
makes the contribution of thermal and suprathermal reaction components 
comparable, while for the D-l-T and D-l-D reactions the effect is rather imper- 
ceptible. One should keep in mind, however, that these results give approx- 
imate picture because the two-temperature model does not reproduce true 
form of particle distributions. Indeed, the knock-on ions are not Maxwellian; 
at least their distribution should be truncated at some critical energy deter- 
mined by kinematics for particle collision. Nevertheless, the above estimations 
indicate that the problem really stands and it is worth studying rigorously. 



3 Plasma kinetic analysis 

In the present work we employ an appropriate plasma kinetic model to study 
the phenomenon on the example of various reactions in self-sustained DT 
fusion plasma with admixture of light elements. The concentration of these 
elements is assumed to be sufficiently low to neglect their role when analyz- 
ing the behavior of main plasma species - fuel ions and 3.5-MeV a-particles 
born in DT fusions. It was shown (see, for example, [1,2,3,4]) that these ener- 
getic a-particles are responsible for non-Maxwellian perturbations of fuel ion 
distributions due to a-D and a-T close collisions. We describe the behavior 
of plasma species a (deuterons, tritons, a-particles) with isotropic velocity 
distributions in terms of a Boltzmann-Fokker-Planck (BFP) equation. The 
BFP equation at steady-state without external heating can be written in the 
following form [4,5]: 



dfc 



dt 



a 



( 



dt 



a 
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C'oul. 
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cond. 



— La + Sa — 



(8) 
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where /„ is the density-normahzed distribution function of species a. The 
plasma is assumed to satisfy the quasi- neutrahty condition: Uf. = n^ + nf + 2nQ. 
The first operator in the right hand of (8) represents the effect of smaU-angle 
a-ion and a-electron Coulomb scattering 




Coul. 



where the functions Aa and Ba are given in [6,7,8]. The second operator is 
a Boltzmann collision integral describing the effect of close a-h coUision (6 = 

0?, a) 



^ Jo Jo 

X / vl^aNEs{v'r)dvl dv'dvl - V fa{v) / Vbfb{Vb) 

\J\v'-v'\ / ^ V Jo 



dt / 

^'^ / NES 



X / v^aNEs{vr)dvr dvb- (10) 

Here v'^ = \v' — v'i,\, Vr = \v — Vb\, P{v' — >• v\vb) gives the probability distribution 
function for the speed v of a scattered particle, and (Tj^es is the collision cross 
section quoted from [9]. The third term 

'dfa\ 1 d ( V-'fa \ ^^^^ 



gives the diffusion in velocity space due to thermal conduction with the typical 
time tIj \ Finally, La and Sa are particle loss and source terms, respectively, 
taking different forms for every ion species. Plasma electrons are considered to 
be Maxwellian at some temperature Tg incorporated in our model by using a 
global power-balance (GPB) equation. This equation determines the relation 
between plasma density and temperature, and has the simple form 

PheatiP'ii Tj) PbremiP'ii '^ei ^^e) Pc.p.iP'ii ^e; ^e) ~ 0. (12) 

Here Pueat is the plasma heating rate by a-particles, Pbrem is the rate of 
bremsstrahlung energy loss, Pep. gives the energy loss due to thermal conduc- 
tion and particle leak. At chosen Tg, being an input parameter in our model, 
plasma density is estimated from (12) assuming 1%^ — na + fit and T^. — 
The detailed description of the kinetic model and explicit expressions for all 
terms in the BFP and GPB equations one can find in [4] . 

Figure 1 shows the particle distributions in energy space calculated under 
conditions close to ITER-like plasma. The non-Maxwellian perturbations of 
the deuteron and triton functions caused by a-D and a-T close collisions are 
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Fig. 1. The energy distribution functions of fuel ions and a-particles calculated 
under two plasma conditions: (i) Te = 10 keV, Tj = 9.8 keV, Ue = 1-2 x 10^'' m~^, 
nd = nt = 5.8 x 10^^ m"^, and (ii) = 20 keV, Tj = 18.8 keV, rie = 6.1 x 10^^ m'^, 
n-d = = 2.8 X 10^^ m~^. The respective Maxwellian distributions are shown by 
the dotted curves. 

clearly marked at energies above a few hundred of keV. ^ The distribution of 
a-particles reflects well slowing down history of these fusion products. This 
distribution reveals moderate energy dependence in the 0.1-3.5 MeV deceler- 
ation range, while at thermal energies it is described by Maxwellian-like form. 
The plasma ion temperature Tj in Fig. 1 is evaluated as 

Ti = (n^"'*=rd + n^""Tt) / (n'^f'' + , (13) 

where the temperatures of deuterons and tritons Tt are obtained by fitting 
the bulk components of the D and T distributions to proper Maxwellian func- 
tions. For conditions considered in the work the difference between and 
does not exceed 4 %. Although it is rather small, we employ the general difini- 
tion of Ti Eq. (13) instead of its reduction Td — Tt — Ti to estimate this basic 
plasma parameter as accurately as possible. An informative parameter is the 



We note that knock-on deuterons were already observed in DT fusion experiments 
at JET [10]. 
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Table 1 

The list of nuclear reactions in the DT /X plasma and the enhancement of their rate 
parameters caused by the knock-on deuterons and tritons 



system 


reaction 


Q-value 




{av)/{av}Mxw 






(MeV) 


(MeV) 


Ti = 10-40 keV 


D+Li 


^Li(d, ni)'Be* 


2.95 


0.429 


2.2 1.5 






4.55 


0.478 


2.5-1.5 




^U{d,pt)a 


2.56 




2.9-1.6 




fix • / 7 \ 

Li(a, a)a 


22.37 




1.5-1.2 


T+Li 




0.51 


0.478 


7.8-1.7 




6Li(t,pi)8Li* 


- 0.18 


0.981 


10^-20 


D+Be 


^Be{d,j)^^B 


15.81 




50-3 


D+C 


i2C(d,pi)i3c* 


- 0.37 


3.089 


10i8_io4 


D+T 


D(t, n)a 


17.59 




< 1.01 


D+D 


D{d, nfRe 


3.27 




< 1.07 



^ Energies of 7 rays emitted by the excited daughter nuclei 

fraction of the knock-on ions n' /n. It is estimated to be at the level of 0.03% 
indicating that the plasma conditions are very close to thermal equilibrium. 



Now we can examine the influence of the knock-on ions on reaction rates in 
the DT/X plasma. In the present study the admixture ions X are chosen to 
be Li, Be, C. These light elements are often considered as low-Z impurity 
in magnetic confinement fusion devices. For example, Li was already used in 
operation of TFTR [11], ASDEX and TEXTOR tokamaks, W-7 AS stellarator 
[12,13], and has been proposed as a diagnostic admixture for fusion machines 
of next generation [14,15,16]. Be and C have been vigorously used for plasma 
diagnostics in JET [17,18]. In order to make the study most informative we 
examine the variety of reactions having different mechanisms. They are hsted 
in Table 1 including some processes proposed for plasma 7-ray spectroscopy, 
the energy-producing Li(ci, a) and tritium-breeding Li{d,pt) reactions, the 
conventional D-|-T and D-|-D fusion processes. The respective reaction cross 
sections are plotted in Fig. 2. 

To calculate (av) for the distorted isotropic distributions of D and T, we 
reduce the general expression for reactivity (2) to the form 

((7v)i2 = / vifiiyi) I V2f2{v2) v^a{v)dvdvidv2, (14) 

711^2 Jo Jo J\vi-V2\ 

where v is the relative speed |vi — V2I. Assuming the comparatively heavy 
particles X to be Maxwellian, the reaction rate parameters have been com- 
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puted at ion temperature Tj = 10-40 keV. The results are plotted in Figs. 3 
and 4, and also displayed in Table 1. We see that although the fraction of the 
knock-on ions is only 0.03%, these particles appreciably affect the D+X and 
T+X reactivity. It is underestimated if Maxwellian DT plasma is assumed, 
and for some reactions the discrepancy between the two approaches becomes 
crucial. The ratio 6 = {crv) / {av)Mxw monotonically increases with decreasing 
Tj down to the plasma ignition point, and for the exothermic reactions pre- 
sented in Fig. 3 it changes approximately within 1.5-2.5 (D-|-Li), 2-8 (T-|-Li) 
and 3-50 (D+Bc). Special attention is worth being paid to the cndothermic 
T-|-Li and D+C reactions displayed in Fig. 4. Here 6 turns out to be several 
orders of magnitude or more, that reflects threshold nature of these processes. 
Both of them are forbidden at energies below thresholds, so only sufficiently 
fast ions contribute to the reactions. At the same time. Fig. 1 shows that the 
amount of these particles is essentially underestimated in Maxwellian plasma. 
The record enhancement is marked for D+C; the high threshold and strong 
Coulomb suppression of the thermal channel make the role of the knock-on 
deuterons extremely important here. Thus, the both reactions proceed via 
suprathermal channels and are solely governed by knock-on ions. This may 
have interesting applications in fusion technology. For example, 0.981-MeV 
photons emitted in ^Li[t,pi) might be applicable to energetic triton and a- 
particle diagnostics [4]. 



In agreement with the comments in Section 1 we find that the knock-on ions 
do not significantly affect the D-l-T and D-l-D reactions. Table 1 shows that 
the enhancement factor 6 does not exceed 1% (D+T) and 7% (D-l-D). 
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Fig. 3. The rate parameters for exothermic reaetions from Table 1. The eomputed 
and MaxwelUan data are shown by the sohd and dotted curves, respectively. The 
dashed curves give the factor S = (crv) / {av)M:x.-w 

4 Conclusion 



We have presented arguments that the Maxwelhan approximation can essen- 
tially underestimate reaction rates in nearly thermal-equilibrium plasma, and 
explicitly demonstrated this phenomenon for some reactions in the DT/X 
plasma. The enhancement of reactivity results from the knock-on perturba- 
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Ion temperature (keV) 



Fig. 4. The rate parameters for endothermic reactions from Table 1. 

tions of ion distributions caused by close collisions with energetic fusion prod- 
ucts. The numerical analysis carried out in the work is consistent with the 
semi-qualitative consideration, but the level of reactivity enhancement marked 
has turned out to be surprisingly high. 

The phenomenon is hkely to have general nature and can play role in other 
plasma systems. This especially concerns threshold nuclear processes - our 
study indicates that evaluation of their rates within the Maxwellian approx- 
imation can involve dramatic errors and prove to be fully useless. It seems 
possible that knock-on ions can affect power balance and even reduce ignition 
temperature for some exotic fuels, in which conditions favorable for the non- 
Maxwellian pumping of ion distributions can be realized. The aneutronic ^He 
plasma would be an interesting object for such a study. Indeed, the ^He-|-^He 
reaction has large Q — 12.9 McV, generates fast charged particles and exhibits 
an appropriate cross section behavior. The reaction cross section rapidly and 
monotonically rises with increasing energy up to the MeV region. 

Apart from laboratory plasmas, the phenomenon discussed may also appear 
in some astrophysical processes. Primordial plasma is of particular importance 
here. Standard big-bang nucleosynthesis (BBN) relies on nuclear reaction net- 
work involving Maxwellian reactivity, and accuracy of nuclear inputs has been 
under attention [19,20]. If the mechanism of non-Maxwellian deviation would 
come into play under BBN specific conditions, it could update input reactivi- 
ties and thereby offer new insight into synthesis of hght elements in the early 
universe. Analysis of this scenario requires coupled cosmological and plasma 
kinetic calculations which were beyond the scope of our work. However, the 
demonstration that a very small, almost negligible fraction of knock-on parti- 
cles (~ 10^"^) in hot plasma can significantly change the rates of some reactions 
gives impetus to such study. 
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